RNA degradome analysis reveals DNE1 endoribonuclease is required for the turnover of diverse mRNA substrates in Arabidopsis.

In plants, cytoplasmic mRNA decay is critical for posttranscriptionally controlling gene expression and for maintaining cellular RNA homeostasis. Arabidopsis DCP1-ASSOCIATED NYN ENDORIBONUCLEASE 1 (DNE1) is a cytoplasmic mRNA decay factor that interacts with proteins involved in mRNA decapping and nonsense-mediated mRNA decay (NMD). There is limited information on the functional role of DNE1 in RNA turnover, and the identities of its endogenous targets are unknown. In this study, we utilized RNA degradome approaches to globally investigate DNE1 substrates. Monophosphorylated 5' ends, produced by DNE1, should accumulate in mutants lacking the cytoplasmic exoribonuclease XRN4, but be absent from DNE1 and XRN4 double mutants. In seedlings, we identified over 200 such transcripts, most of which reflect cleavage within coding regions. While most DNE1 targets were NMD-insensitive, some were upstream ORF (uORF)-containing and NMD-sensitive transcripts, indicating that this endoribonuclease is required for turnover of a diverse set of mRNAs. Transgenic plants expressing DNE1 cDNA with an active-site mutation in the endoribonuclease domain abolished the in planta cleavage of transcripts, demonstrating that DNE1 endoribonuclease activity is required for cleavage. Our work provides key insights into the identity of DNE1 substrates and enhances our understanding of DNE1-mediated mRNA decay.

Integrative Analysis of Gene Expression and miRNAs Reveal Biological Pathways Associated with Bud Paradormancy and Endodormancy in Grapevine.

Transition of grapevine buds from paradormancy to endodormancy is coordinated by changes in gene expression, phytohormones, transcription factors, and other molecular regulators, but the mechanisms involved in transcriptional and post-transcriptional regulation of dormancy stages are not well delineated. To identify potential regulatory targets, an integrative analysis of differential gene expression profiles and their inverse relationships with miRNA abundance was performed in paradormant (long day (LD) 15 h) or endodormant (short day (SD), 13 h) Vitis riparia buds. There were 400 up- and 936 downregulated differentially expressed genes in SD relative to LD budsGene set and gene ontology enrichment analysis indicated that hormone signaling and cell cycling genes were downregulated in SD relative to LD buds. miRNA abundance and inverse expression analyses of miRNA target genes indicated increased abundance of miRNAs that negatively regulate genes involved with cell cycle and meristem development in endodormant buds and miRNAs targeting starch metabolism related genes in paradormant buds. Analysis of interactions between abundant miRNAs and transcription factors identified a network with coinciding regulation of cell cycle and epigenetic regulation related genes in SD buds. This network provides evidence for cross regulation occurring between miRNA and transcription factors both upstream and downstream of MYB3R1.

RNA degradomes reveal substrates and importance for dark and nitrogen stress responses of Arabidopsis XRN4.

XRN4, the plant cytoplasmic homolog of yeast and metazoan XRN1, catalyzes exoribonucleolytic degradation of uncapped mRNAs from the 5' end. Most studies of cytoplasmic XRN substrates have focused on polyadenylated transcripts, although many substrates are likely first deadenylated. Here, we report the global investigation of XRN4 substrates in both polyadenylated and nonpolyadenylated RNA to better understand the impact of the enzyme in Arabidopsis. RNA degradome analysis demonstrated that xrn4 mutants overaccumulate many more decapped deadenylated intermediates than those that are polyadenylated. Among these XRN4 substrates that have 5' ends precisely at cap sites, those associated with photosynthesis, nitrogen responses and auxin responses were enriched. Moreover, xrn4 was found to be defective in the dark stress response and lateral root growth during N resupply, demonstrating that XRN4 is required during both processes. XRN4 also contributes to nonsense-mediated decay (NMD) and xrn4 accumulates 3' fragments of select NMD targets, despite the lack of the metazoan endoribonuclease SMG6 in plants. Beyond demonstrating that XRN4 is a major player in multiple decay pathways, this study identified intriguing molecular impacts of the enzyme, including those that led to new insights about mRNA decay and discovery of functional contributions at the whole-plant level.

Murine Oviductosomes (OVS) microRNA profiling during the estrous cycle: Delivery of OVS-borne microRNAs to sperm where miR-34c-5p localizes at the centrosome.

Oviductosomes (OVS) are nano-sized extracellular vesicles secreted in the oviductal luminal fluid by oviductal epithelial cells and known to be involved in sperm capacitation and fertility. Although they have been shown to transfer encapsulated proteins to sperm, cargo constituents other than proteins have not been identified. Using next-generation sequencing, we demonstrate that OVS are carriers of microRNAs (miRNAs), with 272 detected throughout the estrous cycle. Of the 50 most abundant, 6 (12%) and 2 (4%) were expressed at significantly higher levels (P < 0.05) at metestrus/diestrus and proestrus/estrus. RT-qPCR showed that selected miRNAs are present in oviductal epithelial cells in significantly (P < 0.05) lower abundance than in OVS, indicating selective miRNA packaging. The majority (64%) of the top 25 OVS miRNAs are present in sperm. These miRNAs' potential target list is enriched with transcription factors, transcription regulators, and protein kinases and there are several embryonic developmentally-related genes. Importantly, OVS can deliver to sperm miRNAs, including miR-34c-5p which is essential for the first cleavage and is solely sperm-derived in the zygote. Z-stack of confocal images of sperm co-incubated with OVS loaded with labeled miRNAs showed the intracellular location of the delivered miRNAs. Interestingly, individual miRNAs were predominantly localized in specific head compartments, with miR-34c-5p being highly concentrated at the centrosome where it is known to function. These results, for the first time, demonstrate OVS' ability to contribute to the sperm's miRNA repertoire (an important role for solely sperm-derived zygotic miRNAs) and the physiological relevance of an OVS-borne miRNA that is delivered to sperm.

Analysis of Brachypodium miRNA targets: evidence for diverse control during stress and conservation in bioenergy crops.

BACKGROUND: Since the proposal of Brachypodium distachyon as a model for the grasses, over 500 Bdi-miRNAs have been annotated in miRBase making Brachypodium second in number only to rice. Other monocots, such as switchgrass, are completely absent from the miRBase database. While a significant number of miRNAs have been identified which are highly conserved across plants, little research has been done with respect to the conservation of miRNA targets. Plant responses to abiotic stresses are regulated by diverse pathways many of which involve miRNAs; however, it can be difficult to identify miRNA guided gene regulation when the miRNA is not the primary regulator of the target mRNA. RESULTS: To investigate miRNA target conservation and stress response involvement, a set of PARE (Parallel Analysis of RNA Ends) libraries totaling over two billion reads was constructed and sequenced from Brachypodium, switchgrass, and sorghum representing the first report of RNA degradome data from the latter two species. Analysis of this data provided not only PARE evidence for miRNA guided cleavage of over 7000 predicted target mRNAs in Brachypodium, but also evidence for miRNA guided cleavage of over 1000 homologous transcripts in sorghum and switchgrass. A pipeline was constructed to compare RNA-seq and PARE data made from Brachypodium plants exposed to various abiotic stress conditions. This resulted in the identification of 44 miRNA targets which exhibit stress regulated cleavage. Time course experiments were performed to reveal the relationship between miR393ab, miR169a, miR394ab, and their respective targets throughout the first 36 h of the cold stress response in Brachypodium. CONCLUSIONS: Knowledge gained from this study provides considerable insight into the RNA degradomes and the breadth of miRNA target conservation among these three species. Additionally, associations of a number of miRNAs and target mRNAs with the stress responses have been revealed which could aid in the development of stress tolerant transgenic crops.

Opsin Repertoire and Expression Patterns in Horseshoe Crabs: Evidence from the Genome of Limulus polyphemus (Arthropoda: Chelicerata).

Horseshoe crabs are xiphosuran chelicerates, the sister group to arachnids. As such, they are important for understanding the most recent common ancestor of Euchelicerata and the evolution and diversification of Arthropoda. Limulus polyphemus is the most investigated of the four extant species of horseshoe crabs, and the structure and function of its visual system have long been a major focus of studies critical for understanding the evolution of visual systems in arthropods. Likewise, studies of genes encoding Limulus opsins, the protein component of the visual pigments, are critical for understanding opsin evolution and diversification among chelicerates, where knowledge of opsins is limited, and more broadly among arthropods. In the present study, we sequenced and assembled a high quality nuclear genomic sequence of L. polyphemus and used these data to annotate the full repertoire of Limulus opsins. We conducted a detailed phylogenetic analysis of Limulus opsins, including using gene structure and synteny information to identify relationships among different opsin classes. We used our phylogeny to identify significant genomic events that shaped opsin evolution and therefore the visual system of Limulus We also describe the tissue expression patterns of the 18 opsins identified and show that transcripts encoding a number, including a peropsin, are present throughout the central nervous system. In addition to significantly extending our understanding of photosensitivity in Limulus and providing critical insight into the genomic evolution of horseshoe crab opsins, this work provides a valuable genomic resource for addressing myriad questions related to xiphosuran physiology and arthropod evolution.

The genome of the seagrass Zostera marina reveals angiosperm adaptation to the sea.

Seagrasses colonized the sea on at least three independent occasions to form the basis of one of the most productive and widespread coastal ecosystems on the planet. Here we report the genome of Zostera marina (L.), the first, to our knowledge, marine angiosperm to be fully sequenced. This reveals unique insights into the genomic losses and gains involved in achieving the structural and physiological adaptations required for its marine lifestyle, arguably the most severe habitat shift ever accomplished by flowering plants. Key angiosperm innovations that were lost include the entire repertoire of stomatal genes, genes involved in the synthesis of terpenoids and ethylene signalling, and genes for ultraviolet protection and phytochromes for far-red sensing. Seagrasses have also regained functions enabling them to adjust to full salinity. Their cell walls contain all of the polysaccharides typical of land plants, but also contain polyanionic, low-methylated pectins and sulfated galactans, a feature shared with the cell walls of all macroalgae and that is important for ion homoeostasis, nutrient uptake and O2/CO2 exchange through leaf epidermal cells. The Z. marina genome resource will markedly advance a wide range of functional ecological studies from adaptation of marine ecosystems under climate warming, to unravelling the mechanisms of osmoregulation under high salinities that may further inform our understanding of the evolution of salt tolerance in crop plants.

Heat-induced ribosome pausing triggers mRNA co-translational decay in Arabidopsis thaliana.

The reprogramming of gene expression in heat stress is a key determinant to organism survival. Gene expression is downregulated through translation initiation inhibition and release of free mRNPs that are rapidly degraded or stored. In mammals, heat also triggers 5'-ribosome pausing preferentially on transcripts coding for HSC/HSP70 chaperone targets, but the impact of such phenomenon on mRNA fate remains unknown. Here, we provide evidence that, in Arabidopsis thaliana, heat provokes 5'-ribosome pausing leading to the XRN4-mediated 5'-directed decay of translating mRNAs. We also show that hindering HSC/HSP70 activity at 20°C recapitulates heat effects by inducing ribosome pausing and co-translational mRNA turnover. Strikingly, co-translational decay targets encode proteins with high HSC/HSP70 binding scores and hydrophobic N-termini, two characteristics that were previously observed for transcripts most prone to pausing in animals. This work suggests for the first time that stress-induced variation of translation elongation rate is an evolutionarily conserved process leading to the polysomal degradation of thousands of 'non-aberrant' mRNAs.

AtCCR4a and AtCCR4b are Involved in Determining the Poly(A) Length of Granule-bound starch synthase 1 Transcript and Modulating Sucrose and Starch Metabolism in Arabidopsis thaliana.

Removing the poly(A) tail is the first and rate-limiting step of mRNA degradation and apparently an effective step not only for modulating mRNA stability but also for translation of many eukaryotic transcripts. Carbon catabolite repressor 4 (CCR4) has been identified as a major cytoplasmic deadenylase in Saccharomyces cerevisiae. The Arabidopsis thaliana homologs of the yeast CCR4, AtCCR4a and AtCCR4b, were identified by sequence-based analysis; however, their role and physiological significance in plants remain to be elucidated. In this study, we revealed that AtCCR4a and AtCCR4b are localized to cytoplasmic mRNA processing bodies, which are specific granules consisting of many enzymes involved in mRNA turnover. Double mutants of AtCCR4a and AtCCR4b exhibited tolerance to sucrose application but not to glucose. The levels of sucrose in the seedlings of the atccr4a/4b double mutants were reduced, whereas no difference was observed in glucose levels. Further, amylose levels were slightly but significantly increased in the atccr4a/4b double mutants. Consistent with this observation, we found that the transcript encoding granule-bound starch synthase 1 (GBSS1), which is responsible for amylose synthesis, is accumulated to a higher level in the atccr4a/4b double mutant plants than in the control plants. Moreover, we revealed that GBSS1 has a longer poly(A) tail in the double mutant than in the control plant, suggesting that AtCCR4a and AtCCR4b can influence the poly(A) length of transcripts related to starch metabolism. Our results collectively suggested that AtCCR4a and AtCCR4b are involved in sucrose and starch metabolism in A. thaliana.

Identification of SMG6 cleavage sites and a preferred RNA cleavage motif by global analysis of endogenous NMD targets in human cells.

In metazoans, cleavage by the endoribonuclease SMG6 is often the first degradative event in non-sense-mediated mRNA decay (NMD). However, the exact sites of SMG6 cleavage have yet to be determined for any endogenous targets, and most evidence as to the identity of SMG6 substrates is indirect. Here, we use Parallel Analysis of RNA Ends to specifically identify the 5' termini of decay intermediates whose production is dependent on SMG6 and the universal NMD factor UPF1. In this manner, the SMG6 cleavage sites in hundreds of endogenous NMD targets in human cells have been mapped at high resolution. In addition, a preferred sequence motif spanning most SMG6 cleavage sites has been discovered and validated by mutational analysis. For many SMG6 substrates, depletion of SMG6 resulted in the accumulation of decapped transcripts, an effect indicative of competition between SMG6-dependent and SMG6-independent NMD pathways. These findings provide key insights into the mechanisms by which mRNAs targeted by NMD are degraded.

Differential expression of miRNAs and their target genes in senescing leaves and siliques: insights from deep sequencing of small RNAs and cleaved target RNAs.

MicroRNAs (miRNAs) are a class of small RNAs, which typically function by guiding cleavage of target mRNAs. They are known to play roles in a variety of plant processes including development, responses to environmental stresses and senescence. To identify senescence regulation of miRNAs in Arabidopsis thaliana, eight small RNA libraries were constructed and sequenced at four different stages of development and senescence from both leaves and siliques, resulting in more than 200 million genome-matched sequences. Parallel analysis of RNA ends libraries, which enable the large-scale examination of miRNA-guided cleavage products, were constructed and sequenced, resulting in over 750 million genome-matched sequences. These large datasets led to the identification a new senescence-inducible small RNA locus, as well as new regulation of known miRNAs and their target genes during senescence, many of which have established roles in nutrient responsiveness and cell structural integrity. In keeping with remobilization of nutrients thought to occur during senescence, many miRNAs and targets had opposite expression pattern changes between leaf and silique tissues during the progression of senescence. Taken together, these findings highlight the integral role that miRNAs may play in the remobilization of resources and alteration of cellular structure that is known to occur in senescence.

Sample sequencing of vascular plants demonstrates widespread conservation and divergence of microRNAs.

Small RNAs are pivotal regulators of gene expression that guide transcriptional and post-transcriptional silencing mechanisms in eukaryotes, including plants. Here we report a comprehensive atlas of sRNA and miRNA from 3 species of algae and 31 representative species across vascular plants, including non-model plants. We sequence and quantify sRNAs from 99 different tissues or treatments across species, resulting in a data set of over 132 million distinct sequences. Using miRBase mature sequences as a reference, we identify the miRNA sequences present in these libraries. We apply diverse profiling methods to examine critical sRNA and miRNA features, such as size distribution, tissue-specific regulation and sequence conservation between species, as well as to predict putative new miRNA sequences. We also develop database resources, computational analysis tools and a dedicated website, http://smallrna.udel.edu/. This study provides new insights on plant sRNAs and miRNAs, and a foundation for future studies.

Parallel analysis of RNA ends enhances global investigation of microRNAs and target RNAs of Brachypodium distachyon.

BACKGROUND: The wild grass Brachypodium distachyon has emerged as a model system for temperate grasses and biofuel plants. However, the global analysis of miRNAs, molecules known to be key for eukaryotic gene regulation, has been limited in B. distachyon to studies examining a few samples or that rely on computational predictions. Similarly an in-depth global analysis of miRNA-mediated target cleavage using parallel analysis of RNA ends (PARE) data is lacking in B. distachyon. RESULTS: B. distachyon small RNAs were cloned and deeply sequenced from 17 libraries that represent different tissues and stresses. Using a computational pipeline, we identified 116 miRNAs including not only conserved miRNAs that have not been reported in B. distachyon, but also non-conserved miRNAs that were not found in other plants. To investigate miRNA-mediated cleavage function, four PARE libraries were constructed from key tissues and sequenced to a total depth of approximately 70 million sequences. The roughly 5 million distinct genome-matched sequences that resulted represent an extensive dataset for analyzing small RNA-guided cleavage events. Analysis of the PARE and miRNA data provided experimental evidence for miRNA-mediated cleavage of 264 sites in predicted miRNA targets. In addition, PARE analysis revealed that differentially expressed miRNAs in the same family guide specific target RNA cleavage in a correspondingly tissue-preferential manner. CONCLUSIONS: B. distachyon miRNAs and target RNAs were experimentally identified and analyzed. Knowledge gained from this study should provide insights into the roles of miRNAs and the regulation of their targets in B. distachyon and related plants.

Comprehensive investigation of microRNAs enhanced by analysis of sequence variants, expression patterns, ARGONAUTE loading, and target cleavage.

MicroRNAs (miRNAs) are a class of small RNAs that typically function by guiding the cleavage of target messenger RNAs. They have been shown to play major roles in a variety of plant processes, including development, and responses to pathogens and environmental stresses. To identify new miRNAs and regulation in Arabidopsis (Arabidopsis thaliana), 27 small RNA libraries were constructed and sequenced from various tissues, stresses, and small RNA biogenesis mutants, resulting in 95 million genome-matched sequences. The use of rdr2 to enrich the miRNA population greatly enhanced this analysis and led to the discovery of new miRNAs arising from both known and new precursors, increasing the total number of Arabidopsis miRNAs by about 10%. Parallel Analysis of RNA Ends data provide evidence that the majority guide target cleavage. Many libraries represented novel stress/tissue conditions, such as submergence-stressed flowers, which enabled the identification of new stress regulation of both miRNAs and their targets, all of which were validated in wild-type plants. By combining small RNA expression analysis with ARGONAUTE immunoprecipitation data and global target cleavage data from Parallel Analysis of RNA Ends, a much more complete picture of Arabidopsis miRNAs was obtained. In particular, the discovery of ARGONAUTE loading and target cleavage biases gave important insights into tissue-specific expression patterns, pathogen responses, and the role of sequence variation among closely related miRNA family members that would not be evident without this combinatorial approach.

XRN 5'→3' exoribonucleases: structure, mechanisms and functions.

The XRN family of 5'→3' exoribonucleases is critical for ensuring the fidelity of cellular RNA turnover in eukaryotes. Highly conserved across species, the family is typically represented by one cytoplasmic enzyme (XRN1/PACMAN or XRN4) and one or more nuclear enzymes (XRN2/RAT1 and XRN3). Cytoplasmic and/or nuclear XRNs have proven to be essential in all organisms tested, and deficiencies can have severe developmental phenotypes, demonstrating that XRNs are indispensable in fungi, plants and animals. XRNs degrade diverse RNA substrates during general RNA decay and function in specialized processes integral to RNA metabolism, such as nonsense-mediated decay (NMD), gene silencing, rRNA maturation, and transcription termination. Here, we review current knowledge of XRNs, highlighting recent work of high impact and future potential. One example is the breakthrough in our understanding of how XRN1 processively degrades 5' monophosphorylated RNA, revealed by its crystal structure and mutational analysis. The expanding knowledge of XRN substrates and interacting partners is outlined and the functions of XRNs are interpreted at the organismal level using available mutant phenotypes. Finally, three case studies are discussed in more detail to underscore a few of the most exciting areas of research on XRN function: XRN4 involvement in small RNA-associated processes in plants, the roles of XRN1/PACMAN in Drosophila development, and the function of human XRN2 in nuclear transcriptional quality control. This article is part of a Special Issue entitled: RNA Decay mechanisms.

Methods for validation of miRNA sequence variants and the cleavage of their targets.

MicroRNA (miRNA) variants that share the sequences with other closely related miRNAs have been identified by deep sequencing and have been implicated in the diverse regulation of their target genes. The miRNA variants that originate from the same miRNA precursor are among the most common and have been termed "isomiRs." IsomiRs can be generated by several mechanisms such as differential processing by DICER, RNA degradation, or RNA editing. Members of the same miRNA family that have distinct sequences also contribute to the diversity of miRNA variants. Although many miRNA variants are lowly expressed and may function redundantly with their reference miRNAs, some miRNA variants are highly and/or differentially expressed. In addition, slight differences in sequence among miRNA variants can affect their specificity in target selection. Here, we describe two methods for detecting or validating miRNA variants and the target events they mediate.

Roles of DCL4 and DCL3b in rice phased small RNA biogenesis.

Higher plants have evolved multiple proteins in the RNase III family to produce and regulate different classes of small RNAs with specialized molecular functions. In rice (Oryza sativa), numerous genomic clusters are targeted by one of two microRNAs (miRNAs), miR2118 and miR2275, to produce secondary small interfering RNAs (siRNAs) of either 21 or 24 nucleotides in a phased manner. The biogenesis requirements or the functions of the phased small RNAs are completely unknown. Here we examine the rice Dicer-Like (DCL) family, including OsDCL1, -3a, -3b and -4. By deep sequencing of small RNAs from different tissues of the wild type and osdcl4-1, we revealed that the processing of 21-nucleotide siRNAs, including trans-acting siRNAs (tasiRNA) and over 1000 phased small RNA loci, was largely dependent on OsDCL4. Surprisingly, the processing of 24-nucleotide phased small RNA requires the DCL3 homolog OsDCL3b rather than OsDCL3a, suggesting functional divergence within DCL3 family. RNA ligase-mediated 5' rapid amplification of cDNA ends and parallel analysis of RNA ends (PARE)/degradome analysis confirmed that most of the 21- and 24-nucleotide phased small RNA clusters were initiated from the target sites of miR2118 and miR2275, respectively. Furthermore, the accumulation of the two triggering miRNAs requires OsDCL1 activity. Finally, we show that phased small RNAs are preferentially produced in the male reproductive organs and are likely to be conserved in monocots. Our results revealed significant roles of OsDCL4, OsDCL3b and OsDCL1 in the 21- and 24-nucleotide phased small RNA biogenesis pathway in rice.

MAJOR DEVELOPMENTAL REGULATORS AND THEIR EXPRESSION IN TWO CLOSELY RELATED SPECIES OF PORPHYRA (RHODOPHYTA)(1).

Little is known about the genetic and biochemical mechanisms that underlie red algal development, for example, why the group failed to evolve complex parenchyma and tissue differentiation. Here we examined expressed sequence tag (EST) data from two closely related species, Porphyra umbilicalis (L.) J. Agardh and P. purpurea (Roth) C. Agardh, for conserved developmental regulators known from model eukaryotes, and their expression levels in several developmental stages. Genes for most major developmental families were present, including MADS-box and homeodomain (HD) proteins, SNF2 chromatin-remodelers, and proteins involved in sRNA biogenesis. Some of these genes displayed altered expression correlating with different life history stages or cell types. Notably, two ESTs encoding HD proteins showed eightfold higher expression in the P. purpurea sporophyte (conchocelis) than in the gametophyte (blade), whereas two MADS domain-containing paralogs showed significantly different patterns of expression in the conchocelis and blade respectively. These developmental gene families do not appear to have undergone the kinds of dramatic expansions in copy number found in multicellular land plants and animals, which are important for regulating developmental processes in those groups. Analyses of small RNAs did not validate the presence of miRNAs, but homologs of Argonaute were present. In general, it appears that red algae began with a similar molecular toolkit for directing development as did other multicellular eukaryotes, but probably evolved altered roles for many key proteins, as well as novel mechanisms yet to be discovered.

Massive analysis of rice small RNAs: mechanistic implications of regulated microRNAs and variants for differential target RNA cleavage.

Small RNAs have a variety of important roles in plant development, stress responses, and other processes. They exert their influence by guiding mRNA cleavage, translational repression, and chromatin modification. To identify previously unknown rice (Oryza sativa) microRNAs (miRNAs) and those regulated by environmental stress, 62 small RNA libraries were constructed from rice plants and used for deep sequencing with Illumina technology. The libraries represent several tissues from control plants and plants subjected to different environmental stress treatments. More than 94 million genome-matched reads were obtained, resulting in more than 16 million distinct small RNA sequences. This allowed an evaluation of ~400 annotated miRNAs with current criteria and the finding that among these, ~150 had small interfering RNA-like characteristics. Seventy-six new miRNAs were found, and miRNAs regulated in response to water stress, nutrient stress, or temperature stress were identified. Among the new examples of miRNA regulation were members of the same miRNA family that were differentially regulated in different organs and had distinct sequences Some of these distinct family members result in differential target cleavage and provide new insight about how an agriculturally important rice phenotype could be regulated in the panicle. This high-resolution analysis of rice miRNAs should be relevant to plant miRNAs in general, particularly in the Poaceae.

MicroRNAs as master regulators of the plant NB-LRR defense gene family via the production of phased, trans-acting siRNAs.

Legumes and many nonleguminous plants enter symbiotic interactions with microbes, and it is poorly understood how host plants respond to promote beneficial, symbiotic microbial interactions while suppressing those that are deleterious or pathogenic. Trans-acting siRNAs (tasiRNAs) negatively regulate target transcripts and are characterized by siRNAs spaced in 21-nucleotide (nt) "phased" intervals, a pattern formed by DICER-LIKE 4 (DCL4) processing. A search for phased siRNAs (phasiRNAs) found at least 114 Medicago loci, the majority of which were defense-related NB-LRR-encoding genes. We identified three highly abundant 22-nt microRNA (miRNA) families that target conserved domains in these NB-LRRs and trigger the production of trans-acting siRNAs. High levels of small RNAs were matched to >60% of all ∼540 encoded Medicago NB-LRRs; in the potato, a model for mycorrhizal interactions, phasiRNAs were also produced from NB-LRRs. DCL2 and SGS3 transcripts were also cleaved by these 22-nt miRNAs, generating phasiRNAs, suggesting synchronization between silencing and pathogen defense pathways. In addition, a new example of apparent "two-hit" phasiRNA processing was identified. Our data reveal complex tasiRNA-based regulation of NB-LRRs that potentially evolved to facilitate symbiotic interactions and demonstrate miRNAs as master regulators of a large gene family via the targeting of highly conserved, protein-coding motifs, a new paradigm for miRNA function.